Scaling function for the noisy Burgers equation in the soliton approximation 
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We derive the scaling function for the one dimensional noisy Burgers equation in the two-soliton 
approximation within the weak noise canonical phase space approach. The result is in agreement 
with an earlier heuristic expression and exhibits the correct scaling properties. The calculation 
presents the first step in a many body treatment of the correlations in the Burgers equation. 



The strong coupling aspects of driven noncquilibrium 
systems present an important challenge in statistical 
physics. The phenomena in question are widespread, in- 
cluding turbulence, interface and growth problems, and 
chemical reactions. 

Here the noisy Burgers equation or equivalently the 
Kardar-Parisi-Zhang (KPZ) equation, describing the 
growth of an interface, is one of the simplest models of a 
driven system showing scaling and pattern formation. 

In one dimension the noisy Burgers equation for the 
slope u = Vh of a growing interface has the form |Qj 

— - AuV \u = vV 2 u + V?y , (1) 
at J 

(r){xt)ri(x't')) = AS{x - x')6(t - t') . (2) 

The height h is then governed by the KPZ equation [|| 

dh/dt = vV 2 h + (\/2)(Vh) 2 + n . (3) 

In (|l|) and (^) v is the damping, A the nonlinear mode 
coupling, and r\ a Gaussian white noise of strength A, 
correlated according to (0). The equation (|l|) is more- 
over invariant under the Galilean transformation 



Xuot 



U + Uq 



(4) 



The Burgers equation ([!]) and its KPZ equivalent in one 
and higher dimensions and related models in the same 
universality class have been studied intensely in recent 
years owing to their importance as models for a class of 
nonequilibrium systems ^ |j . 

We have in a series of papers analyzed the one dimen- 
sional case defined by ([!]) and (Q) in an attempt to un- 
cover the physical mechanisms underlying the pattern 
formation and scaling behavior fj]. Emphasizing that 
the noise strength A is the relevant nonperturbative pa- 
rameter, driving the system into a stationary state, the 
method was initially based on a weak noise saddle point 



approximation to the Martin-Siggia-Rose functional for- 
mulation of the noisy Burgers equation. This work 
was a continuation of earlier work based on the mapping 
of a solid-on-solid model onto a continuum spin model 
||. More recently the functional approach has been su- 
perseded by a canonical phase space method || deriving 
from the symplectic structure |HJ of the Fokker-Planck 
equation associated with the Burgers equation. 

The functional or the equivalent phase space approach 
valid in the weak noise limit A — > yields coupled deter- 
ministic mean field equations 



- AuV ) u= vV 2 u - V 2 p , 

dt 



— - AuV ) p= -vVp 



(5) 
(6) 



for the slope u and a canonically conjugate noise field 
p (replacing the stochastic noise 77), determining orbits 
in a canonical phase space and replacing the stochastic 
Burgers equation (Q). The equations @ and @) derive 
from a principle of least action with Hamiltonian density 



H = p{vV 2 u+\u\Ju- 
with an orbit v! — > W 



S(u",t, 



-(l/2)V 2 p) and action S associated 
traversed in time t, 



t,u 



, du 

I V-Qt 



H 



(7) 



The action is of central importance and serves as a 
weight for the nonequilibrium configurations (cp. the 
Boltzmann-Gibbs factor exp(—E/kT) for equilibrium 
systems). The action moreover gives access to the time 
dependent and stationary probability distributions 



P{u , t, u) oc exp [— S{u , t, u )/A] , 



(8) 



and P st {u") = lim^oo P(u", t, u') , and associated mo- 
ments, e.g., the slope correlations 

(uu){xt) = [ J[du u"(x)u'(0)P(u",t,u')P st (u') . (9) 
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The equations (||) and (||) admit static soliton solutions 

= /iutanh[k s x] , k s = \uj2v , fi = ±1 , (10) 

moving solitons are generated by the Galilean boost (|J). 
Denoting the right and left boundary values by u+ and 
U-, respectively, the velocity is given by 



= -2v/X 



(11) 



The index \i labels the right hand soliton for n = 1 on 
the 'noiseless' manifold p = 0, also a solution of the 
damped noiseless Burgers equation for r\ = 0; and the 
noise-excited left hand soliton for fi = — 1 on the 'noisy' 
manifold p — 2vu, a solution of the growing (unstable) 
noiseless Burgers equation for v — > —v. The wavenumber 
k s sets the inverse length scale. The field equations also 
admit linear mode solutions superimposed on the soli- 
tons; for A = they become the usual diffusive modes of 
the driven equation ||. 

The physical picture emerging from this analysis is a 
many body formulation of the pattern formation of a 
growing interface in terms of a dilute gas of propagating 
solitons with superimposed linear modes. The formula- 
tion also associates energy E = J dxH and momentum 
IT = J dxuVp with a soliton mode, yielding the disper- 
sion law 



E oc (A/^ 1/2 )IT 



(12) 



with dynamic exponent z = 3/2 and it follows that the 
strong coupling fixed point features are associated with 
the soliton dynamics, i.e., defect or domain wall excita- 
tions. 

In this Letter we pursue the form of the slope corre- 
lations @; the basic building block in the many body 
formulation. We focus in particular on the scaling func- 
tion which is of central importance. The dynamic scal- 
ing hypothesis |||| and general arguments based on the 
renormalization group fixed point structure (^Jll[l imply 
the following long time-large distance form of the slope 
correlations in the stationary state: 



(uu)(xt) = (A/2i/)x- 2(1 " c) G(a;/e(i)) 



(13) 



Here G is the scaling function and £ = 1/2 the roughness 
exponent inferred from the known stationary probability 
distribution 



P st (w) oc exp 



-(f/A) / dx u 



(14) 



Within the canonical phase space approach ( pT| ) follows 
from the structure of the zero-energy manifold which at- 
tracts the phase space orbits for t — ► oo. The dynami- 
cal exponent z = 3/2 then follows from the scaling law 
£ + z = 2 implied by the Galilean invariance (^) [^). 
In the present approach the exponent z is inferred from 



the (gapless) soliton dispersion law (|12[). Finally, the 
growth of lateral correlations along the interface is char- 
acterized by the time dependent correlation length £(i). 
In the nonlinear nonequilibrium Burgers case de- 
scribes the propagation of solitons and is given by £(i) = 
(A/is) 1 ^ 3 (At) 2 / 3 . In the linear equilibrium Edwards- 
Wilkinson case |t3| £(t) characterizes the growth of diffu- 
sive modes and has the form £(t) = (vt) 1 / 2 . It moreover 
follows from the 'fluctuation-dissipation theorem' (|l4| ) 
that u is uncorrelated and that the static correlations 
have the form, (uu)(x) = (A/2v)6(x), independent of A. 
This is consistent with £ = 1/2 and the limiting form of 
the scaling function lim^—^ G(w) = 1 for x 3> £(t). In 
the dynamical regime for £(t) 3> x the correlation decay, 
i.e., (uu)(x,t) — > (u)(u) — 0, and the scaling function 
vanishes like G(w) oc u' 2 ^ 1 "''-' for w — > 0. 

Generally, the scaling function can be inferred from 
the asymptotic properties of the slope correlations (|^). 
In order to evaluate those we must i) determine an or- 
bit from u' to u" in time t by solving the field equa- 
tions (j^) and (^|) as an initial-final value problem in u 
(p is a slaved variable), ii) evaluate the associated action 
(0) in order to weigh the orbit and determine P (note 
that P s t is given by (|l4|)), and, finally, iii) integrate over 
initial and final configurations v! and u". Even in the 
one dimensional case discussed here such a calculation 
appears rather formidable in the general multi-soliton - 
linear mode case and we must resort to partial results. 

In the weak noise limit the action ^ according to (||) 
provides a selection criterion determining the dominant 
dynamical configuration contributing to the distribution. 
For A — > an important contribution to the growth mor- 
phology is constituted by two-soliton configurations or 
pair excitations 

U2(x, t) — u+(x — vt — xi) + u~(x — vt — x 2 ) , (15) 

obtained by matching two Galilei-boosted static solitons 
of opposite parity (// = ±1) ( ^0| ) centered at x\ and x-i 
with soliton separation t = \x\ — xi \ and amplitude 2u. 
According to the soliton condition ([ll]) the pair excita- 
tion propagates with velocity v = — Xu and has vanishing 
slope field u = at the boundaries, corresponding to a 
horizontal interface. 

Whereas the solitons ( |l0| ) for /i = ±1 lie on the 
transient and stationary submanifolds (separatrixes) for 
p s = (the 'noiseless' kink) and p s = 2vu s (the 'noisy' 
kink), respectively, and constitute the 'quarks' in the 
many body formulation, the pair excitation (^), satisfy- 
ing the boundary conditions, is the elementary excitation 
or 'quasi particle' (in the Landau sense) in the present 
scheme and is characterized by the dispersion law (|T^). 
The pair excitation is an approximate solution to the field 
equations (g) and (^) with a finite lifetime |l4|]. Over a 
time scale controlled by the damping v the pair decays 



2 



into diffusive modes; this is consistent with the observa- 
tion that the phase space orbits approach the zero-energy 
manifold for t — > oo. 

Unlike a general multi-soliton configuration which 
changes in time owing to soliton-soliton collisions, the 
pair excitation preserves its shape over a finite time pe- 
riod, see ref. P|l4|]. Imposing periodic boundary condi- 
tions for the slope field the motion of a pair with ampli- 
tude 2u corresponds to a simple growth mode where the 
height field h, i.e., the integrated slope field, increases 
layer by layer for each revolution of the soliton pair in a 
system of size L. From the KPZ equation (||) it follows 
that (dh/dt) — (X/2)(u 2 ) in a stationary state. Setting 
u — » 2u this is consistent with the increase Ah = 2u£ 
during the passage time At — t/v — t/\u for a soliton 
pair of size £. In Fig. [j] we have depicted the two-soliton 
growth mode in the slope field u and the associated height 
field h. The pair excitation, which can also be conceived 
as a bound state composed of two solitons, has the am- 
plitude 2u, size £, carries energy E 2 = — (16/3)^ A|m| 3 , 
momentum IT2 = — 4z/u|w|, and action 



S 2 = -vX\u\ 3 t 



(16) 



Using the definition ([)]) it is an easy task to evaluate 
the contribution to the slope correlations from a single 
pair. The normalized stationary distribution P st is ob- 
tained from ( pr]) by insertion of (|l5|). Considering the 
inviscid limit for v — > we have 



P at {u,l) = fVt 1 exp [-4(v/A)u 2 £] 
O st = {it A/ v) 1 1 2 1? 1 2 . 



(17) 
(18) 



Correspondingly, inserting (|l6|) in (j|) the normalized soli- 
ton pair transition probability is 

P sol (u, t) = fi- o { cxp [-(4/3)(v/A)X\u\ 3 t\ , (19) 
Q sol = (2/3)r(l/3)[(3/4)(A/ i ,)(l/At)] 1 /3 . (20) 

We note that the normalization factor f2 st for the sta- 
tionary distribution varies as L 3 ! 2 and that the distri- 
bution thus vanishes in the infinite size limit; moreover, 
the mean size of a pair is equal to L, characteristic of an 
extended excitation (a string). Likewise the transition 
probability P so \ goes to zero for large times in accordance 
with the decay of a soliton pair into diffusive modes. 

The evolution of (uu)(xt) in the two-soliton sector is 
straightforward. The final configuration u" is simply the 
initial configuration u' displaced vt along the axis with no 
change of shape, i. e., u"{x) = u'(x + vt), v = —Xu'. Not- 
ing that the integral over u' and u" — u' only contributes 
when the pair configurations overlap and integrating over 
the size £ we obtain the slope correlations 

£ Q Jdue-^n e -^^B Cl (u) 
(uu){xt) = — _ 4 r— — , ( 21 ) 



Jdue-^-hc 2 {i 



where the cut-off functions originating from the overlap 
are given by C\(u) = l/Au 2 - (1 + l/4u 2 ) exp(— 4u 2 ) and 
C 2 (w) = (l/4w 2 )(l — exp(— 4w 2 )), respectively. In order 
to facilitate the discussion of (21) we have introduced 
the noise- induced length and time scales £0 — A/v and 
to = A/vX; note that A = ^oAch and, moreover, the 
crossover or saturation time t s — to(L/£ ) 3 / 2 ; the cor- 
relation length is then £ = £o(t/to) 2 / 3 . The expression 
([H]) holds for t > and is even in x (seen by changing u 
to —it). It samples the soliton pair propagating with ve- 
locity Ait and is in general agreement with spectral form 
discussed in the 'quantum' treatment in ||. In Fig. || we 
have shown the two-soliton overlap configurations con- 
tributing to the slope correlations. 

The weight of single soliton pair is of order 1/L and 
the correlation function (uu) thus vanishes in the ther- 
modynamic limit L — > 00. For a finite system L enters 
setting a length scale together with the saturation time 
t s cx L 3 / 2 defining a time scale, and (uu) is a function 
of x/L and t/t s as is the case for the two-soliton expres- 
sion (^l|). This dependence should be compared with 
the wavenumber decomposition of (uu) for A = 0. Here 
(uu)(xt) oc (l/L) Ylnjto ex P(~ (2im) 2 t/ L 2 ) exp(innx/L), 
depending on x/L and t/L 2 , corresponding to the satu- 
ration time t s oc L 2 , z = 2. Keeping only one mode for 
n = 1 (uu) has the same structure as in the soliton case. 
In the linear case we can, of course, sum over the totality 
of modes and in the thermodynamic limit L — > 00 replace 
(1/L) J2 n by J dk/2n obtaining the intensive correlations 
(uu)(xt) = (A/2i/)(4tt^)- 1 / 2 exp(— x /2vt). Similarly, 
we expect the inclusion of multi-soliton modes to allow 
the thermodynamic limit to be carried out yielding an 
intensive correlation function in the Burgers case. 

For a finite system we have in general |l5| (uu) (xt) = 
(l/L)G L {x/L,t/L 3 / 2 ) with scaling limits: G L (x/L,0) cx 
const, for x ~ L, Gl(x/L, 0) oc L/x for x L and 
G L (0,t/L 3 / 2 ) oc const, for t > L 3 / 2 , G L (0, i/L 3 / 2 ) oc 
Lit 2 ' 3 for t <C L 3 / 2 . For L — > 00 we obtain 
G L (x/L,t/L 3 ' 2 ) -> {L/x)G(x/t 2 / 3 ) in conformity with 
©■ 

It is an important feature of the two-soliton expression 
( pl| ) that the dynamical soliton interpretation directly 
implies the correct dependence on the scaling variables 
x/L and t/t s oc i/L 3 / 2 , i.e., independent of a renor- 
malization group argument. However, the scaling lim- 
its are at variance with Gl- Setting, according to (|2l| ) 
(uu)(xt) = (£ /L)F(x/ L,t/t s ), F(x/L,0) assumes the 
value .47 for x <C L and decreases monotonically to the 
value ~ .08 for x ~ L, whereas Gl diverges as L/x 
for x <§; L. Likewise, F(0,t/t s ) decays from .47 for 
t < t s oc L 3/2 to for t > t s ; for t ~ t s we have F ~ .15, 
whereas Gl diverges as L/t 2 / 3 for t < t s . 

This discrepancy from the scaling limits is a feature the 
two-soliton contribution which only samples the correla- 
tion from a single soliton pair. Moreover, at long times 
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the soliton contribution vanishes and the scaling function 
is determined by the diffusive mode contribution in ac- 
cordance with the convergence of the phase space orbits 
to the stationary zero-energy manifold. We note, how- 
ever, the general trend towards a divergence for small 
values of x and t is a feature of F. 

Introducing the scaling variables w — x/£ oc x/t 2 / 3 
and t == t/t s oc tjl?! 2 we can also express ( f2l| ) in the 
form 



(uu){xt) = {£ /L)F 2 (w,t) , 
where the scaling function F is now given by 



(22) 



F(w, t) = ± Tr-,3 — , 23 

and summarize our findings in Fig. || where we have de- 
picted F(w, t) for a range of r values. For fixed small 
w = x/£ oc xjt 1 ! 3 we have F — > .47 for r = f/f s oc 
t/L 3 / 2 -> 0; for large r we obtain i* 1 — > 0. The motion of 
the weak maximum towards smaller values of w for de- 



creasing r is a feature of the functional form of F in (23), 
i.e., the soliton approximation, and probably not a prop- 
erty of the true scaling function which is not expected to 
have any particularly distinct features ||-|{| . 

In this Letter we have presented the two-soliton con- 
tribution to the slope correlations and ensuing scaling 
function within the weak noise canonical phase space ap- 
proach to the noisy Burgers equation in one dimension. 
The expression is in accordance with a general spectral 
form proposed earlier on the basis of the many body in- 
terpretation of a growing interface and has the correct 
scaling dependence. This calculation presents the first 
step in a many body or field theoretical treatment of the 
correlations in the noisy Burgers equation based on a 
transparent physical quasi particle picture of the growth 
mechanisms and ensuing morphology. Details will be pre- 
sented elsewhere. 

Discussions with A. Svane, J. Hertz, B. Derrida, M. 
Lassig, J. Krug and G. Schiitz are gratefully acknowl- 
edged. 
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